Magnetic fluctuations from stripes in cuprates 
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Within the time-dependent Gutzwiller approximation for the Hubbard model we com- 
pute the magnetic fluctuations of vertical metallic stripes with parameters appropriate for 
Lai.875Bao.i25Cu04 (LBCO). For bond- and site-centered stripes the excitation spectra are sim- 
ilar, consisting of a low-energy incommensurate acoustic branch which merges into a "resonance 
peak" at the antiferromagnetic wave vector and several high-energy optical branches. The acoustic 
branch is similar to the result of theories assuming localized spins whereas the optical branches are 
significantly different. Results are in good agreement with a recent inelastic neutron study of LBCO. 

PACS numbers: 74.25.Ha,71.28.+d,74.72.-h,71.45.1r 
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Since the discovery of cuprate superconductors the elu- 
cidation of their magnetic properties has been a subject 
of intense research in the high-T c community due to their 
possible relevance for the superconducting mechanism |lj. 
However, no consensus has been reached yet wether the 
magnetic excitation spectra in the various cuprate mate- 
rials can be traced back to some universal phenomenol- 
ogy which could be expected in the face of the robust 
nature of superconductivity. The insulating parent com- 
pounds show long-range antiferromagnetic (AF)order in 
the Cu02 planes below the Neel temperature J2|. This 
static AF order is lost above a concentration of added 
holes per planar copper x ~ 0.02 but complex dynamical 
spin correlations persist up to the overdoped regime [lllil- 

Regarding the lanthanum cuprates (LCO), neutron 
scattering (NS) experiments have revealed low-energy 
incommensurate magnetic excitations at wave-vectors 
(hereafter in units of 2n/a 4J) (1/2 ± e, 1/2) and 
(1/2, 1/2 ± e)[3 in the doping regime where the material 
is superconducting (x > 0.05). The incommensurabil- 
ity e depends linearly on doping e = x up to x ~ 1/8 
and stays constant beyondl3|. This behavior can be 
successfully explained^ 0, 13 by metallic, self-organized 
quasi-one-dimensional structures called stripes[9j which 
are oriented parallel to the Cu-0 bond (hereafter vertical- 
stripes) and act as antiphase domain walls for the anti- 
ferromagnetic order. The strongest experimental support 
for the stripe scenario stems from neutron|l0Ull| and x- 
ray scattering [12j experiments on LCO where co-dopants 
induce a pinning of the stripes so that the associated 
charge order can be detected. 

In YBa 2 Cu30 6+y (YBCO) incommensurate magnetic 
fluctuations have been detected by inelastic NS (INS) 13, 
[lj, |la, 116j with a similar doping dependence of the low- 
energy incommensurability to that of LCO[l5|. Upon 
increasing energy the incommensurate branches contin- 
uously disperse towards the so-called resonance mode 
at wave-vector Qaf = (1/2, 1/2) a collective magnetic 
mode that grows up below T c [l3, lla, ll2|- The energy 



E r of the spin resonance seems to scale linearly with 
T c which has led to speculations that it could be re- 
lated to the superconducting pairing and phase coher- 
ence (cf. Ref. J2fJ,|2l| and references therein), thus sug- 
gesting a magnetic origin of the high transition tem- 
peratures of cuprates. Finally above E r the magnetic 
fluctuations in YBCO acquire again an incommensurate 
structure PuL UM. l22l 123 . |24| and are already observed 
above T c [lfJ. 

Very recent experiments were dedicated to explore 
the problem of universality in the magnetic excita- 
tions between different cuprate families J23J, 12.4 I26L 
l27j . In particular Tranquada and collaborators |27J re- 
ported INS measurements of the magnetic excitations in 
Lai.875Barj.i25Cu04 which shows static charge and spin 
order. Most interestingly the dispersion of spin exci- 
tations shows features which resemble closely those of 
YBCO although the measurement has been performed 
above T c . The similarity between both compounds has 
been further demonstrated by a high energy study of 
YBCO 16] and a high resolution experiment on opti- 
mally doped LCOpfif. These works[3 IJJ lil Efil Wfy 
are certainly an important step towards a unified under- 
standing of magnetic fluctuation in cuprates. On the 
other hand they open new questions. The spin exci- 
tations in YBCO are commonly explained in terms of 
an itinerant picture as arising from a dispersing two- 
particle bound state induced by AF correlations in a 
d-wave superconducting system|25j, |28(. On the other 
hand magnetic excitations on top of stripes are usually 
described in a localized moment picture within spin- 
wave theory (LSWT)@ H3- If the origin of mag- 
netic excitations in cuprates is universal, what picture 
is more appropriate? The spectra of LSWT are difficult 
to reconcile with several features of YBCO [la, 123, l25| 
and LC Ol26l |27| challenging the stripe interpretation 
itself |23j, |26J. Here we show that the magnetic excita- 
tions in Lai.875Bao.i25Cu04 (and possibly in the other 
cuprates) can be understood in terms of the spin flue- 



tuations of metallic stripes which are in an intermediate 
regime, i.e. neither the localized nor the itinerant picture 
applies. Contrary to the LSWT computations |29j our re- 
sults are in agreement with experiment 27] over the whole 
range of energies and momenta. 

Our investigations are based on the one-band Hubbard 
model with hopping restricted to nearest (~ t) and next 
nearest (~ t') neighbors 
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Here c\ a destroys (creates) an electron with spin a at site 

i, and n, )CT = c\ a Ci y(T . U is the on-site Hubbard repulsion. 

Static charge and spin textures are obtained within an 
unrestricted Gutzwiller approximation (GA). Dynamical 
properties are computed on top of the inhomogeneous 
solutions within the time-dependent GA|3lJ (TDGA) 
which previously has been shown to provide an accu- 
rate description of the optical conductivity within the 
more realistic three-band model|32j. Here, because we 
restrict to low-energy magnetic excitations, we do not 
expect that interband effects will play an important role 
and therefore a one-band description should be sufficient. 

Parameters were fixed by requiring that a) the lin- 
ear concentration of added holes is l/(2a) according to 
experiment [2I [lfj, |2JJ and b) a TDGA computation of 
the undoped AF insulator reproduces the experimental 
dispersion relation|33J. Condition a) was shown to be 
very sensitive to t'/t [8j whereas condition b) is sen- 
sitive to U/t and t, the former parameter determining 
the observed energy splitting between magnons at wave- 
vectors (1/2,0) and (1/4, 1/4) 0. Indeed, the splitting 
vanishes within spin-wave theory applied to the Heisen- 
berg model which corresponds to U/t — > oo. We find 
that both conditions are met by t'/t = —0.2, U/t = 8 
and t = 354meV. In this way all parameters are fixed 
and the subsequent computation of the magnetic fluctu- 
ation of the stripes can be considered without free pa- 
rameters. Results shown below are obtained in a 40 x 40 
sites lattice and for doping x = 1/8 corresponding to a 
period of charge modulation d = 4a. The charge and 
spin structure of the static GA solution is similar to the 
one reported in Ref. |8j where a slightly different param- 
eter set has been used. The local magnetization of BC 
stripes is reduced to 0.20 (0.38) with respect to 1/2 in- 
side (outside) the core of the charge domains. The band 
structure is quasi-one dimensional with a Fermi momen- 
tum Uf ~ 1/8. Bond-centered (BC) and site-centered 
(SC) stripes are almost degenerate in energy which is 
also found in other one-band calculations joL [34 and 
within the three-band model at higher doping[7|,|32j. We 
anticipate that for both structures the magnetic excita- 
tion spectra are rather similar [see Fig. dc)] with differ- 
ences regarding the intensity distribution and gaps be- 
tween magnetic bands. For definiteness we mainly re- 
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FIG. 1: 10 contour levels of ujS ± (uj,q) at integer values show- 
ing the dispersion of magnetic excitations along (a) and per- 
pendicular (b) to an array of BC stripes. The factor of 
u> serves to cancel an infrared intensity divergence. Panel 
(c) shows S(u), q) integrated on the magnetic Brillouin zone 
around (1/2,1/2) 37] for BC and SC stripes together with the 
experimental data 27] . For the meaning of "error" bars see 
Ref. \gft. 

strict ourself to the BC case because of marginally bet- 
ter agreement with experiment (see below). Additionally 
these textures constitute the more stable configuration at 
x = 1/8 in the more accurate three-band model 7] and in 
first principle computations j35j . 

We compute the transverse dynamical structure factor, 
S -L( w ,q) = ^ B ^ -£ v \(0\S+\u)\ 2 6(w-u v ) which is 
probed by INS. Here g — 2 and as in LSWT we included 
the renormalization factor Z x 2] . In the insulating phase 
our spectral weights are close to LSWT and therefore 
we adopt Z x — 0.51 133|. r\ takes into account polariza- 
tion factors in the NS cross section and will be discussed 
below. Energy is sampled at intervals of 3.5meV. 

A contour level plot of w^fw, q) [Fig.[IJa),(b)] shows 
the dispersion of magnetic excitations for stripes ori- 
ented along the y-axis. The lower (acoustic) branch 
perpendicular to the stripes (b) is similar to the lowest 
branch found within LSWT[29|. Indeed, since the acous- 
tic branch involves long-wavelength excitations it should 
not depend on the short-range details of the model. It 
shows the correct Goldstone-like behavior going to zero 
frequency at the ordering wave- vector Q s = (1/2 ±e, 1/2) 
with e = 1/8. Starting from Q s one observes two 
branches of spin-waves where the one dispersing towards 
smaller q x rapidly loses intensity. The other one remains 
very intense up to Qaf where it can be associated with 
the resonance peak. Moreover, as shown in Fig. QJa,b) 
the dispersion develops a local maximum at Qaf ex- 
plaining the strong intensity in the momentum integrated 
structure factor [Fig. E^c)] a * the resonance frequency. 
In the direction of the stripe the excitations display a 
"roton-like" minimum. The energy of the resonance at 
Qaf is E r = 65meV which is reasonably close to the 



experimental one for this system E r = 55meV 27] taking 
into account the simplicity of the model and the fact that 
we have no free parameters. 

In order to compare with experiments one should av- 
erage over the two possible orientations of the stripes. 
In Fig. ^ this amounts to adding to each of the panels 
(a) and (b) the data of the other panel reflected with 
respect to the central axes. One obtains an X-shaped 
dispersion providing a natural explanation for the X-likc 
feature seen in both YBCO[Il 0,111 and LCOH^. 

The acoustic branch and its continuation in Fig. Q^a) 
are quite similar to the dispersion obtained in a weakly 
coupled two-leg ladder system |38| . In these approaches 
parameters are adjusted to drive the system into the 
quantum critical point separating the quantum param- 
agnet from the magnetically ordered state |3SJ. Hence 
the similarity with our spectra is reasonable since re- 
garding the magnetism the ordered state corresponds to 
the present ground state and one can expect continuity 
of the excitations at the transition. The spin-leg lad- 
der theories 38], however, rely heavily on the fine tuning 
of the coupling parameter between the legs and an even 
charge periodicity of the stripes (like for e = 1/8). In con- 
trast we obtain qualitatively similar spectra for e = 1/8 
and e = 1/10 (x — 0.1, d = 5a) 0, |g|, in accord with 
experiment 26] . Note that this does not exclude the in- 
teresting possibility that a small spin-gap opens by the 
ladder mechanism for x > 1/8 and e = 1 /8 |35l I38L I39J . 

Above the acoustic branch in Fig. ^b) there are three 
optical branches. The lowest two almost touch at q x — 
1/4. Gap positions are in agreement with LSWT but the 
dispersion is not. Indeed the computations in Ref. |29j 
yield a shift of the two lower optical branches by 1/8 
along q x with respect to our computation. Formally this 
difference can be traced back to the non local nature of 
magnetic excitations in our system which shows magnetic 
moments considerably smaller than in the insulator, and 
more importantly is metallic. In LSWT only processes 
of the kind Sf S~ are allowed in the effective interaction 
kernel with i and j being close neighbors 29]). We find 
that not only processes with i far from j are important 
but also processes of the kind e] i-cy^ej i cy^f with all four 
sites different. 

The absolute spectral weights are determined by 
the parameter 77. If one assumes a three-dimensional 
isotropic distribution of magnetic domains 77 = 2/3. If 
instead the magnetization is parallel to the CuO plane 
77 will decrease with energy and take values closer to 
1/2. We adopt for 77 80% of the isotropic value which 
gives good account of the lower part of the spectra. (At 
higher energies a further decrease is expected [23 )■ For 
a slightly more doped sample |36| we find 77 ~ 2/3. 

Some details of the relative spectral weight are better 
reproduced by BC stripes, i) The resonance energy E r 
is closer to the experimental value in the BC case, ii) 
For BC stripes the gap between the acoustic branch and 



the first optical branch produces a dip in the integrated 
intensity at 80meV followed by a peak at lOOmeV which 
is in agreement with experiment contrary to the SC case. 
Although these features favor BC stripes more theoretical 
and experimental work is needed to confirm this finding. 

Fig. reports constant frequency scans of 5 , - L (w, q). 
At the lowest energy for y-axis oriented stripes [panel (a) 
left] one intersects the Goldstone mode close to the in- 
commensurate wave- vector Q s . At higher energy [(b), (c)] 
one intersects the spin-wave cone. The intensity is very 
anisotropic along the intersection and only the region 
closest to Qaf is visible. This effect is also apparent from 
Fig. Q] Although our study focusses on an underdoped 
system the large intensity difference fits nicely with the 
fact that only the branch closest to Qaf is seen in a re- 
cent high resolution study on optimally doped LCO[26j. 
Immediately below the resonance [Fig. El (d)] the spin 
excitations have shifted close to Qaf and one observes 
the appearance of intensity in the stripe direction due to 
the "roton-like" behavior. At 81meV one intersects the 
first gap of Fig. ^ (b) so that the main features are due 
to excitations along the stripes. At higher energies the 
contributions from the modes perpendicular to the stripe 
become gradually more important. Moreover, the orien- 
tation averaged scans are in remarkably good agreement 
with the corresponding panels of Ref. [27J • 

Another very interesting finding are the horizontal 
cigar like features seen at 35meV and higher energies. 
Similar features have been obtained in weak coupling 
in the longitudinal channel 40] . These structures ap- 
pear close to q y ~ 1/4 ~ 2kp and correspond to spin- 
flip backscattering processes of the quasi-one dimensional 
metallic subsystem. At this q y the continuum extends to 
about lOOmeV. As expected there is also a sharp collec- 
tive mode corresponding to forward scattering processes 
but with too weak intensity to be observable in Fig. [2] 
Diffusive scattering revealed a similar phenomenon in 
La5/3Ni!/ 3 04 where stripes are insulating but have a 
spin- 1/2 degree of freedom at the core|4l|. 

Concluding, we have calculated the magnetic excita- 
tions of cuprates in the stripe phase within the TDGA 
applied to an extended Hubbard model. From a method- 
ological point of view the present work together with our 
previous study Ref. [32| shows that TDGA is accurate 
and simple enough in order to obtain realistic dynam- 
ical properties of textured strongly correlated system. 
Stripes in our computation are metallic charge and spin 
density- waves in a regime intermediate between localized 
moment and itinerant pictures. Our results are in good 
agreement with experiments in Lai.875Bao.i25Cu04 pro- 
viding a straightforward explanation of the energy and 
momentum dependent evolution of S ± (q,oj) in terms of 
stripes and hence confirming stripes as a robust feature of 
this system with a firm theoretical basis. The fact that 
the magnetic excitation spectra are similar in different 
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FIG. 2: Constant frequency scans of S x (w,q) for BC stripes convolved with a Gaussian (FWHM=0.06 r.l.u.). The title 
indicates the energy or the energy window (e-i) over which intensities are averaged. For each figure (a)-(i) the left panel shows 
the intensity for stripes oriented along the y-direction (divided by 2). In order to compare with the experimental data of 
Ref. |23 the right panel shows the average of both horizontal and vertical stripes with the coordinate system rotated by 45° 
with wave-vector units of 2-7r/(\/2a). Parameters are as in Fig.^and as explained in the text. 



cuprate families and for different doping levels suggests 
that stripes or the proximity to stripe instabilities are a 
universal property of cuprates and hence that they may 
be relevant for the superconducting mechanism. 
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